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(54) Method and apparatus to minimize phase misregistration artifacts in gated CT images 



(57) A method and apparatus for generating an op- 
timally registered working set of two-dimensional CT im- 
ages in a gated CT system are described. The method 
includes, where an acquisition period is divided into se- 
quential acquisition phases, obtaining CT image data for 
separate phase ranges during each acquisition phase, 
dividing the ROI into several adjacent sub-volumes gen- 



erating phase image sets including 2D images for each 
phase range within each sub-volume and selecting one 
phase range image set within each sub-volume corre- 
sponding to each acquisition phase to generate a work- 
ing image set where the selected phase range image 
sets are optimally aligned. 
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Description 

[0001 J The field of the invention is gated computer- 
ized tomography (CT) imaging and more specifically 
methods and apparatus that reduce the amount of mis- 
alignment between parallel two dimensional images in 
a working CT image set. 

[0002] Many different types of medical imaging sys- 
tems have been developed that are used for different 
purposes. Perhaps the most common type of imaging 
system category includes X-ray systems wherein radi- 
ation is directed across a portion of a patient to be im- 
aged and toward a detector panel. An exemplary X-ray 
detector panel includes a Csl scintillator coupled with an 
amorphous silicon array. With radiation directed toward 
a region of a patient to be images (i.e., a region of inter- 
est), the region of interest blocks some of the radiation 
and some of the radiation passes through the region and 
is collected by the panel. The amount of radiation that 
passes through the region along the trajectory of a given 
radiation ray depends upon the type of tissue along the 
trajectory. Thus, a tumor may block more radiation than 
flesh and bone may block more radiation than a tumor 
and so on. Hence X-ray system can be used to collect 
a "projection" through a patient. 

[0003] While useful, simple X-ray systems have many 
limitations. One important limitation to X-ray imaging 
systems is that such systems, as described above, only 
provide side projections through a region and cannot be 
used to generate other useful images such as "slice" im- 
ages (i.e., images perpendicular to projection images) 
through a region of interest. For instance, an exemplary 
useful slice image may include a slice image through a 
patient's heart. 

[0004] Another type of imaging system that is useful 
in generating slice images is generally referred to as a 
computerized tomography (CT) system. An exemplary 
CT system includes a radiation source and a radiation 
detector mounted on opposite sides of an imaging area 
where the imaging area is centered along a translation 
or Z-axis. The source generates radiation that is colli- 
mated into a beam including a plurality of radiation rays 
directed along trajectories generally across the imaging 
area. A line detector may be positioned perpendicular 
to the Z-axis to collect slice image data during a data 
acquisition period. 

[0005] During an acquisition period a region of inter- 
est is positioned within the imaging area and, with the 
radiation source turned on, the region of interest blocks 
some of the radiation and some of the radiation passes 
through the region and is collected by the line detector. 
As in X-ray systems, the amount of radiation that passes 
through the region of interest along the trajectory of a 
given radiation ray is dependent upon the type of tissue 
along the trajectory. In CT systems the source and line 
detector are rotated about the region of interest within 
a rotation plane through the region of interest so that 
radiation "projections" can be collected for a large 



number of angles about the region. By combining the 
projections corresponding to a slice through the region 
of interest using a filtering and back projecting tech- 
nique, a two-dimensional tomographic or axial image (i. 

5 e., a slice image) of the slice is generated. 

[0006] While some diagnostic techniques only require 
one or a small numberof slice images, many techniques 
require a large number of parallel CT slice images. For 
example, some techniques require examination of many 

10 parallel images to identify where an arterial blockage be- 
gins and ends and the nature of the blockage there be- 
tween. As another example, many techniques reformat 
two dimensional data into, in effect, three dimensional 
volumetric images, that can be sliced and diced in sev- 
eral different directions so that various image planes can 
be employed. For instance, where two dimensional data 
is acquired for transverse or cross sectional slices 
through a three dimensional region of interest (e.g., 
through a patient's torso), the data may be reformatted 
20 to generate sagital (i.e., the side plane passing through 
the long axis of the body) or coronal (i.e., the frontal 
plane passing through the long axis of the body) images 
through the region of interest. 

[0007] In order to generate several slice images rap- 

25 idly, CT detectors are typically configured having sever- 
al parallel detector rows such that, during a single rota- 
tion about the imaging area, each detector row collects 
data that can subsequently be used to generate a sep- 
arate CT slice image. 

30 [0008] While increasing the number of detector rows 
reduces acquisition time, detector elements are relative- 
ly expensive and thus more rows translates into a more 
costly overall system. As a balance between cost and 
speed, most multi-row detectors include less than 1 0 de- 

35 tector rows. Hereinafter it will be assumed that an ex- 
emplary detector includes eight detector rows. 
[0009] Where a detector includes eight rows and 
more than eight slice images are required, several dif- 
ferent acquisition periods are typically used to acquire 

40 the necessary slice image data. For instance, assume 
that 80 slice images (an admittedly small number but 
sufficient for exemplary purposes) through a ROI are re- 
quired. In this case, the ROI may be divided into ten sep- 
arate sub-volumes, each of the ten sub-volumes corre- 

45 sponding to a separate eight of the 80 required slice im- 
ages. Thereafter, ten separate acquisition periods may 
be used to collect the sets of slice image data corre- 
sponding to the ten sub-volumes, data corresponding to 
eight separate slice images collected during each of the 

50 ten separate acquisition periods. 

[0010] It has been found that, for large volumes or 
ROIs, data necessary to generate many parallel thin 
slice images can be acquired most rapidly by helically 
collecting the data. To this end, while the source and 

55 detector are rotated about the imaging area, a patient 
bed is translated there through so that the radiation fan 
beam sweeps a helical path through the ROI. After hel- 
ical data is collected, the data is converted to slice image 
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data by any of several different weighting and filtering 
processes and thereafter the slice image data is back- 
projected to form the viewable image. 
[001 1 ] In the case of helically acquired and stored raw 
data, the data can be used to construct virtually any 
number of slice images through a corresponding ROI. 
For instance, despite using a detector having eight rows 
of elements to collect helical data, the data may be proc- 
essed to generate 1 6 ; 20, 500 or even thousands of sep- 
arate slice images or : indeed, may be interpolated to 
generate a 3-D volumetric image, if desired. 
[0012] In most imaging systems that generate still im- 
ages, it is important that, to the extent possible, during 
data acquisition, the structure being imaged remain 
completely still. Even slight structure movement during 
acquisition can cause image artifacts in, and substan- 
tially reduce the diagnostic value of, resulting images. 
For this reason, during acquisition periods, patients are 
typically instructed to maintain the region of interest 
within the imaging area as still as possible by, for in- 
stance, holding the patient's breadth. 
[0013] Despite a patient's attempts to control move- 
ment, certain anatomical structures cannot be held still 
and continue movement during acquisition periods. For 
instance, a patient's heart beats continually during data 
acquisition cycles and the beating movement compli- 
cates the process of acquiring diagnostic quality data. 
[001 4] In the case of the heart, fortunately, the beating 
cycle is repetitive and there are certain cycle phases 
during which the heart muscle is relatively at rest. As 
well known in the art, during a diastolic phase of the 
beating cycle when the heart is filling with blood, the 
heart is relatively at rest and movement is minimal. 
Thus, by restricting data acquisition periods to the di- 
astolic phases of the heart beating cycle, relatively 
movement-free data can be acquired and used to gen- 
erate CT slice images. 

[001 5] To this end, the industry has developed cardiac 
gated CT imaging systems. These systems generally 
take two different forms including shoot and move gating 
scans and retro-gating reconstructions. In thexase of 
shoot and move scans, an electrocardiogram (EKG) 
system is used to monitor heart beating phase and to 
gate the acquisition of data so that data is only acquired 
during specific phases of the heart beating cycle (e.g., 
systolic, diastolic, etc.) Thereafter the data is used to 
generate slice images in a conventional manner. In the 
case of retro-gating reconstruction, a full set of helical 
data is acquired and stored along with corresponding 
EKG signals. Thereafter, a heart cycle phase range is 
selected which indicates a range of the cycle for which 
images should be generated and an image reconstruc- 
tor retrieves the helical data sub-set corresponding to 
the phase range from each heart cycle and generates 
the required images. 

[001 6] In addition to minimizing movement related im- 
age artifacts, each of the gating processes (i.e., pro- 
spective and retrospective) is also meant to reduce mis- 



registration between sets of images that are generated 
using data corresponding to different sub-volumes of a 
region of interest. For instance, in the case above where 
a region is divided into ten separate sub-volumes and 
5 data for each sub-volume is collected during a separate 
acquisition period, if data for two consecutive sub-vol- 
umes were collected during different heart beating 
phase, resulting images would likely be misaligned. 
Thus, by collecting data for all sub-volumes during sim- 
10 iiar heart beating phases, misalignment is substantially 
reduced. In the cases of axially acquired data and heli- 
cally acquired data this means restricting data to a spec- 
ified phase range within each heart beating cycle. For 
instance, the acquired period may be between 70% and 
'5 80% of the total heart beating cycle where the cycle be- 
gins and ends at peak cycle amplitudes. 
[0017] Hereinafter the phrase "phase location" will be 
used to refer to as a phase point within a heart cycle and 
the phrase "phase range" will be used to refer to a range 
20 that is centered on a corresponding phase location. 
[0018] Unfortunately, despite cardiac gating efforts, it 
has been recognized that misregistration or misalign- 
ment of sub-volume images can still occur for several 
reasons. First, as known in the industry, EKG signals 
25 provide only an indirect way to measure heart motion 
and therefore cannot be used to precisely identify iden- 
tical phase locations within a heart beating cycle. Sec- 
ond, it is known that, while generally periodic, the heart 
muscle may not go through precisely the same motions 
30 during consecutive heart cycles so that, even if precise 
phase locations could be identified within a heart beat- 
ing cycle, those locations may not correspond to a sim- 
ilarly positioned heart. Third, in the case of high heart 
rates (i.e., a child's heart) the gating system may have 
35 inadequate temporal resolution to facilitate proper gat- 
ing. These gating problems are further exacerbated 
when attempting to generate images including coronary 
arteries as segments of a given artery may be at rest at 
somewhat different phase locations within consecutive 
40 heart beating cycles. 

[0019] Gating related phase mis-registrations can be 
quite apparent in ventricle walls when viewing an image 
dataset corresponding to a multi-planar reformat ren- 
dering from a sagital or coronal perspective. Similarly, 
45 the mis-registrations are apparent in the coronary arter- 
ies when viewing an image dataset with a curved refor- 
mat rendering. 

[0020] It has been recognized that when gated CT 
techniques are employed to collect image data corre- 

50 sponding to adjacent ROI sub-volumes during sequen- 
tial acquisition phases, where each acquisition phase is 
further sub-divided into shorter phase ranges, often im- 
age sets corresponding to the different phase ranges 
within sequential acquisition phase align better than im- 

55 age sets corresponding to the same phase ranges with- 
in sequential acquisition phases. For example, assum- 
ing adjacent first and second sub-volumes correspond- 
ing to first and second sets of eight slice images and 
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that data corresponding to the first and second sub-vol- 
umes is to be collected during the diastolic phases of 
first and second heart beating cycles, respectively. In 
this case, each of the first and second diastolic phases 
may be divided into beginning and ending phase ranges 
including first and second halves of the first and second 
diastolic phases /respectively. 

[0021] Also, assume that during the first diastolic 
phase two sets of image data for the first sub-volume (i. 
e. ( for the eight slices of the first sub-volume) are ob- 
tained, a first set obtained during the beginning phase 
range and a second set obtained during the ending 
phase range. Similarly, assume that during the second 
diastolic phase two sets of image data for the second 
sub-volume (i.e., forthe seven slices of the second sub- 
volume) are obtained, a first set obtained during the be- 
ginning phase range and a second set obtained during 
the ending phase range. Even though the beginning 
phase range data sets correspond to similar phase 
ranges of the heart beating cycle : it has been found that 
images generated using the beginning phase range da- 
ta sets are often characterized by greater mis-registra- 
tion than images generated using the beginning phase 
range data set from the first diastolic phase and the end- 
ing phase range data set from the second diastolic 
phase. 

[0022] Thus, it has been recognized that where data 
is collected for each sub-volume of a ROI during sepa- 
rate diastolic phases, instead of obtaining a single data 
set corresponding to each diastolic phase, the diastolic 
phase can be divided into several phase ranges and a 
separate sub-volume data set can be obtained for each 
phase range and thereafter, during post acquisition 
processing, image sets corresponding to the different 
phase ranges can be compared and the sets that align 
or register most accurately can be combined into a work- 
ing image set for further diagnostic purposes. While the 
example above where each diastolic phase is divided 
into beginning and ending phase ranges facilitates bet- 
ter results than systems that do not divide the diastolic 
phase, alignment is generally further enhanced as the 
number of diastolic phase divisions are increased. For 
instance, generally, dividing each diastolic phase into 
five phases ranges yields better results than dividing 
each diastolic phase Into two phase ranges. 
[0023] Throughout this specification, while some of 
the examples are described in the context of either a 
prospective gating method or a retrospective gating 
method, it should be understood that the present inven- 
tion is useable with in prospective or retrospective meth- 
ods and processors and should not be limited to one or 
the other. It should suffice to say that where an example 
is taught in the context of one or the other types of sys- 
tems, the type of system not literally taught has only 
been omitted in the interest of simplifying this explana- 
tion and not to limit the invention in any way. 
[0024] It should also be noted that the present inven- 
tion is also useful in the case of advanced multi-sector 



reconstruction algorithms to improve temporal imaging 
resolution. These algorithms and how the present inven- 
tion would be used therewith should be obvious to one 
of ordinary skill in the art in light of the specification 
5 which follows. 

[0025] The invention will now be described in greater 
detail, by way of example, with reference to the draw- 
ings, in which:- 

*o Fig. 1 is a perspective view of a CT apparatus used 
to practice the present invention which includes a 
detector array having rows and columns of detector 
elements and fan beam source; 

Fig. 2 is a block diagram of CT control system which 
may be used to control the CT apparatus of Fig. 1 
and which is useful for the purposes of practicing 
the present invention; 

Fig. 3 is a schematic diagram illustrating a region of 
interest subdivided into sub-volumes and a two-di- 
mensional image stack corresponding thereto; 

Fig. 4 is a graph and schematic diagram illustrating 
an exemplary heart beaing cycle, an exemplary di- 
astolic phase and phase ranges corresponding to 
the diastolic phase and also image sets corre- 
sponding to each of the phase ranges; 

Fig. 5 is similar to Fig. 4 albeit illustrating a series 
of diastolic phases and corresponding image sets; 

Fig. 6 is a schematic diagram illustrating two-dimen- 
sional images in coronal view which correspond to 
different diastolic phases and generally the same 
phase range within each of the diastolic phases; 

Fig. 7 is a schematic diagram illustrating a portion 
of the diagram in Fig. 6; and 

Fig. 8 is a flow chart illustrating a method according 
to the present invention. 

A. Hardware 

[0026] Referring now to Fig. 1 , a CT scanner for use 
with the present invention includes a gantry 20 having 
an opening that defines an imaging area (not separately 
numbered) where gantry 20 supports an x-ray source 
10 oriented to project a fan beam 40 of x-rays along a 
beam axis 41 through a patient 42 to an opposed detec- 
tor array 44. The gantry 20 rotates to swing the beam 
axis 41 within a gantry plane 38 defining the x-y plane 
of a Cartesian coordinate system. Rotation of gantry 20 
is measured by beam angle B from an arbitrary refer- 
ence position within the gantry plane 38. 
[0027] A patient 42 resets on a patient support table 
46 which may be moved along a translation axis 48 
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aligned with a Z-axis of the Cartesian coordinate sys- 
tem. Table 46 crosses gantry plane 38 and is radio- 
translucent so as not to interfere with the imaging proc- 
ess. 

[0028] The x-rays of the fan beam 40 diverge from the 5 
beam axis 41 within the gantry plane 38 across a trans- 
verse axis 50 generally orthogonal to both the beam axis 
41 and the translation axis 48 at a fan beam angle y. 
The x-rays of beam 40 also diverge slightly from the 
beam axis 41 and the gantry plane 38 across the trans- 
lation axis 48 (i.e., along the Z axis). 
[0029] After passing through patient 42, the x-rays of 
the fan beam 40 are received by detector array 44 which 
has detector elements 18' arranged in eight rows ex- 
tending along the traverse axis 50 and a plurality of col- 
umns extending along the translation axis 48. The sur- 
face of detector array 44 may be planar or may follow a 
section of a sphere or cylinder having a center at focal 
spot 26 or alternatively at the system isocenter. 
[0030] The detector elements 1 8' each receive x-rays 
and provide intensity measurements along separate 
rays of the fan beam 40. Each intensity measurement 
describes the attenuation via a line integral of one fan 
beam ray passing through a portion of volume or region 
of interest (ROI) 43 of patient 42. ROI 43 is typically wid- 
er along the Z-axis than the slice volume measured by 
a conventional CT system fan beam along the Z-axis. 
The rows of detector elements 1 8* subdivide the detec- 
tor array and hence the fan beam along the Z-axis. 
[0031] Referring now to Fig. 2, an exemplary control 
system for controlling the CT imaging system of Fig. 1 
includes a plurality of gantry associated control modules^ 
52, a table motor control 58, a computer 60, an opera- 
tor's console 65 and a mass storage device 66. The gan- 
try associated control modules 52 include an x-ray con- 
trol module 54.. a gantry motorcontrol module 56, a data 
acquisition system 62 and an image reconstructor 68. 
The x-ray control 54 provides power and timing signals 
to the x-ray source 10 to turn the source on and off as 
required under the control of computer 60. The gantry 
motor control 56 controls the rotational speed and posi- 
tion of the gantry 20 and provides information to com- 
puter 60 regarding gantry position. The table motor con- 
trol 58 controls translation speed of table 46 and pro- 
vides position feedback information back to computer 
60. 

[0032] Data acquisition system 62 samples and digi- 
tizes intensity signals from the detector elements 1 8' of 
detector array 44 and provides the digitized signals to 
computer 60 which in turn stores the signals in mass 
storage device 66. A slip ring connects all gantry mount- 
ed elements to other system components that are not 
mounted tc the gantry for two way communication as 
well known in the art. After data is collected, image re- 
constructor 52 is controlled to combine the collected da- 
ta to form images for display via console 65 or some 
other display device. 

[0033] Referring still to Figs. 1 and 2, computer 60 



runs a pulse sequencing program to perform the inven- 
tive data processing method as described in more detail 
below. To this end, computer 60 receives commands 
and scanning parameters via operator console 65 which 
is generally a CRT display and keyboard. Console 65 
allows an operator to enter parameters for controlling a 
data acquiring scan, to select images to be displayed 
and to display reconstructed image and other informa- 
tion from computer 60. A mass storage device or mem- 
ory 66 provides a means for storing operating programs 
for the CT imaging system, as well as image data for 
future reference by the operator. Both computer 60 and 
image reconstructor 52 have associated electronic 
memory (not shown) for storing data. 
[0034] In operation, gantry motor control 56 brings 
gantry 20 up to a rotational speed and table motorcon- 
trol 58 begins translation of the table 46. The x-ray con- 
trol 54 turns on the x-ray source 10 and projection data 
is acquired on the continuous basis as the table is 
moved through, and the gantry 20 is rotated about, the 
imaging area. At each beam angle B, the projection data 
acquired comprises intensity signals corresponding to 
each detector element 18' at each particular column and 
row of array 44. After data acquisition, the data is stored 
as helical data in mass storage device 66 and can be 
weighted and filtered to generate slice image data cor- 
responding to separate trans-axial slice images through 
the ROI 43 using any of various processes as well 
known in the art. 

B. Assumptions 

[0035] Referring now to Figs. 1 and 4, for the purpos- 
es of this explanation, it will be assumed that ROI 43 
includes patient 42's heart characterized by a heart 
beating waveform 60 that has a typical cycle period P 
between peak amplitude points R1 and R2 and that has 
a diastolic phase DP that generally occurs between 
heart beating phase times t-j and x 8 . . In addition, it will 
be assumed that helically acquired data is stored in its 
raw helical format initially and is only converted to slice 
image data after a system operator selects a heart beat 
phase at which to view an image of the heart. 
[0036] In addition, while the invention is applicable to 
either retrospective or prospective gating systems, the 
invention will be described herein in the context of a pro- 
spective shoot and move gating system where a proc- 
essor controls data acquisition so that acquisition is only 
acquired during diastolic periods. 
[0037] Moreover, referring also to Fig. 3, it will be as- 
sumed that ROI 43 has a dimension Lr along the Z-axis 
48 and that a typical helical scanning rate only enables 
collection of data corresponding to a z-axis dimension 
Ls during a diastolic acquisition phase DP where dimen- 
sion Ls is one tenth of the dimension Lr. Thus, ROI 43 
is dividable into ten separate sub-volumes 122, 124, 
126, 128... 130, 132 and 134 (only seven shown), a sep- 
arate sub-volume corresponding to each diastolic 
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phase. 

[0038] Furthermore, unless indicated otherwise, it will 
be assumed that a system operator requests that the 
system generate a total of 70 equi-spaced transaxial 
slice images through ROI 43. To this end, in Fig. 3, an 
exemplary slice image stack is identified by numeral 1 04 
with separate stack end images identified by numerals 
1 00 and 1 02. Thus, with ten separate sub-volumes and 
a required 70 equi-spaced slice images, seven separate 
images correspond to each sub-volume. In Fig. 3 the 
separate image sets corresponding to sub-volumes 
122, 124, 126, 128, 130, 132 and 134 are identified by 
numerals 108, 108, 110, 112, 114, 116 and 118, respec- 
tively. In Fig. 4, the diastole acquisition phase DP illus- 
trated corresponds to sub-volume 128 and image set 
112 generally. 

[0039] In addition, it will be assumed that, while the 
diastolic phase ranges from times t 1 through t 8 (see Fig. 
4), enough data is collected during a shorter phase 
range to generate the images within a sub-volume im- 
age set. For instance, referring again to Figs. 3 and 4, 
during a phase range v N between times t 3 and ^ data 
sufficient to generate slice images in set 1 1 2 is acquired. 

C. Theory 

[0040] Referring now to Fig. 5, a plurality of sequential 
heart beating cycles C1 , C2,. . C1 0, like the cycle in Fig. 
4 are illustrated, where each cycle includes a separate 
diastolic phase DP1 , DP2,... DP10, respectively. In ad- 
dition, referring to Figs. 3 and 5, separate image sets 
106, 108... 118 corresponding to sub-volumes 122, 
124... 1 34 are also illustrated. Consistent with prior EKG 
gating techniques, during each diastolic phase, data 
corresponding to a separate set of images and a sepa- 
rate sub-volume is collected. For instance, during first 
phase DP1 data corresponding to transaxial slice image 
set 106 (i.e., including images 1-7) is collected, during 
phase DP2 data corresponding to transaxial slice image 
set 108 (i.e., including images 8-14) is collected and so 
on. 

[0041] Referring still to Fig. 4, after the helical data is 
stored and during a post-acquisition period, when 
choosing images to be viewed, a system operator se- 
lects a specific phase point or location SP from within 
the diastolic phase DP at which the image reconstructor 
(see 68 in Fig. 1 ) should create the 70 required transax- 
ial slice images. Thereafter, reconstructor 68 retrieves 
a sub-set of image data corresponding to phase ranges 
centered on the phase location SP within each diastolic 
phase and uses the data sub-sets to generate the 70 
required images. For instance, in Fig. 4, with selected 
phase location SP as illustrated, reconstructor 68 se- 
lects the data sub-set corresponding to phase range v N 
to generate slice image set 1 06(v N ). Similarly, referring 
to Fig. 5, reconstructor 68 selects the data sub-set cor- 
responding to phase range v N within diastolic phase 
DP2 to generate a slice image set 108(v N ) (not sepa- 



rately illustrated), selects the data sub-set correspond- 
ing to phase range v N within diastolic phase DP10 to 
generate a slice image set 1 1 8(v N ) (not separately illus- 
trated) and so on. Thereafter reconstructor 68 employs 
5 any of several different well known processes to gener- 
ate slice images corresponding to the separate selected 
data sub-sets. 

[0042] Referring now to Figs 1 and. 6, a schematic 
illustrates exemplary image sets within the X-Y transax- 

w ial plane that may be generated by reconstructor 68 in 
the manner described above using data corresponding 
to phase ranges v N . Consistent with the description 
above, image sets are identified by numerals 106v N , 
108v N , etc., with each set corresponding to a separate 

15 ROI sub-volume (e.g., 122, 124, etc. in Fig. 3). While 
seven separate slice images correspond to each image 
set, separate images are only illustrated for image set 
• 112 (v N ). 

[0043] In Fig. 6, the outline of an anatomical structure 

20 151 (e.g., the left ventricle wall of patient 42's heart) is 
captured by the image sets. It can be seen that some 
adjacent image sets are aligned such that the structure 
outline 151 appears continuous between adjacent im- 
age sets. For instance, as illustrated, image set 106v N 

25 appears to be aligned with image sets 1 08v N and 1 1 0v N 
and image set 116v N appears to be aligned with image 
sets 114v N and 108v N . However, because of structure 
movement and other system operating nuances, struc- 
ture outline segments corresponding to other adjacent 

so image sets are misaligned or misregistered. For in- 
stance, one misregistration occurs between image sets 
110v N and 112v N . Other misregistrations are illustrated 
in Fig. 6 and the misregistrations are collectively identi- 
fied by numeral 200. 

35 [0044] When the slice images illustrated in Fig. 6 are 
further manipulated to generate sagital or coronal imag- 
es or even thick axial images that combine slice images 
from different sub-volumes, the misregistrations 200 
cause artifacts that reduce diagnostic value appreciably. 

40 [0045] Ithasbeen recognized that, in many cases, im- 
age sets corresponding to different phase ranges within 
sequential diastolic phases align better than image sets 
corresponding to the same phase ranges within sequen- 
tial diastolic phases. To this end, referring to Fig. 4, di- 

45 astolic acquisition period DP may be sub-divided into 
separate phase ranges v N _ 2 , v H _ v v N , v N+1 and v N+2 
where the duration of each phase range corresponds to 
sufficient helical data to generate a corresponding slice 
image set. For instance, data to generate seven image 

so slice data set 112(v N+2 ) corresponds to phase range 
v N+2» data to generate seven image slice data set 112 
(v N+1 ) corresponds to phase range v N+1 , data to gener- 
ate seven image slice data set 112(v N _.,) corresponds to 
phase range v N ^ and data to generate seven image 

55 slice data set 112(v N _ 2 ) corresponds to phase range 
v n-2 "Thus, in the present example, 35 separate images 
arranged in five separate phase dependent sets 112 
(v N - 2 ). 112(v N . 1 ),112(v N ),112(v N+1 ) and 112(v N+2 ) can 
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be generated during diastolic phase DP where each set 
corresponds to the same sub-volume (e.g., 128 in the 
present example as illustrated in Fig. 3). 
[0046] Referring also to Fig. 7, a schematic corre- 
sponding to section 171 of Fig. 6 is illustrated and in- 
cludes the left portions of structure outline 1 51 that cor- 
respond to image sets 1 1 0v N and 1 1 2v N and segments 
of the structure outline corresponding to image sets 
108v N and 113v N . In addition, section 171 also illus- 
trates the left portions of structure outline 151 that cor- 
respond to image sets 112(v N+1 ) and 112(v N _ 2 ) (the left 
portions of outline 1 51 corresponding to image sets 1 1 2 
(v N+2 ) and 112(v N . 1 ) are not illustrated in Fig. 7). Clearly 
image set 112(y N . 2 ) registers (i.e., is aligned) more 
closely with adjacent image set 110(v N ) than does im- 
age set 112(v N ) and therefore misregistration between 
the image sets generally can be reduced if, assuming 
that image set 110(v N ) is included in a working image 
set, image set 112(v N _ 2 ) is substituted for image set 112 
(v N ). 

[0047] Similarly, image set 112(v N+1 ) registers (i.e., is 
aligned) precisely with adjacent image set 110(v N ) and 
therefore misregistration between the image sets gen- 
erally can be reduced if, again assuming image set 1 1 0 
(v N ) is included in a working image set, image set 112 
( v n+i) is substituted for image set 112(v N ). 
[0048] Thus, by generating a plurality of phase range 
dependent image sets for each diastolic phase and then 
identifying best aligned image sets, an optimal working 
image set can be generated that has better diagnostic 
characteristics. In the present example, referring again 
to Fig. 4, during each diastolic phase DP1, DP2, etc., 
the period is divided into five separate phase ranges 
v n-2> v n-i> v n» v n+i and v n+2 and *h e helical phase range 
dependent image data is then used to generate five sep- 
arate image sets where each set includes seven sepa- 
rate images. Thereafter, the image sets corresponding 
to adjacent ROI sub-volumes (see again 128 - 134 in 
Fig. 3) are compared to identify an optimal working set 
including a separate image set corresponding to each 
sub-volume where the optima! set may include image 
sets from different phase ranges. 

D. Operation 

[0049] Referring now to Fig. 8, a retrospective gating 
method 140 according to the present invention is illus- 
trated. Referring also to Figs. 1, 2 and 3, at process 
block 1 69, a system operator uses operator console 65 
to identify a ROI 43 within patient 42 through which 70 
transaxial two-dimensional slice images 104 should be 
generated. 

[0050] At block 144, during data acquisition, the im- 
aging system 38 collects helical CT image data for the 
entire ROI 43 and stores the helical data as a function 
of EKG signals (i.e., as a function of acquisition cycles) 
within mass storage device 66. 

[0051] After all of the helical data has been acquired 



for ROI 43 and stored in storage device 66, control pass- 
es to block 146 where a system operator used console 
65 to identify a phase command that indicates a partic- 
ular phase location to be imaged. To this end, referring 

5 again to Fig. 4, the exemplary selected phase location 
is identified as SP. Continuing, at block 148, computer 
60 identifies separate phase ranges from within each 
diastolic phase as a function of the phase location SP. 
In the present example, referring still to Fig. 4 and also 

io to Fig. 5, computer 60 sub-divides each diastolic phase 
into five separate phase ranges where a central range 
v N is centered on selected location SP and is between 
times x 3 and x 6 two leading phase ranges and two fol- 
lowing phase ranges that lead and follow the centered 

'5 phase range v N , respectively. The leading phase ranges 
include a phase range <j> N _ 2 ranging between times x 1 
through x 4 and phase range v N-1 that extends between 
times x 2 and t 5 . The following phase ranges include 
phase range v N+1 ranging between times x 4 through t 7 

20 and phase range (|> N+2 that occurs between times x 5 and 
x 8 . Thus, in the present example, the phase ranges 
overlap somewhat. 

[0052] Continuing, at block 142, computer 60 re- 
trieves data from device 66 corresponding to each of the 

25 phase ranges and for each of the diastolic phases (i.e., 
retrospective gating function) and at block 1 52 computer 
60 causes image reconstructor 68 to generate phase 
image sets including 2D images for each phase range 
for each diastolic period. For instance, referring once 

30 again to Fig. 4, image, reconstructor 68 generates a sep- 
arate image set including seven images for each of 
phase period v N _ 2 , v N _.,, v N , v N+1 and v N+2 for sub-vol- 
ume 128 (see also Fig. 3) where sub-volume 128 cor- 
responds to one of the diastolic phases. Similarly, refer- 

35 ring also to Fig. 5, five separate phase specific image 
sets are generated for each of the other sub-volumes 
122, 124, 126, ... 130, 132 and 1 34 where each of these 
sub-volumes corresponds to a separate one of the di- 
astolic phases. All of the image sets are stored within 

40 mass storage device 66. 

[0053] Continuing, and again referring to Figs. 2, 4, 5 
and 6, computer 60 selects the phase image sets that 
are most closely aligned with the selected phase loca- 
tion SP and generates a working image set including all 

45 of the selected phase image sets. In the present exam- 
ple, the phase image set most closely aligned with lo- 
cation SP correspond to phase range v N and, as de- 
scribed above are illustrated in Fig. 6. This step occurs 
at block 1 54. At block 156, computer 60 compares all of 

so the working set image sets (see again Fig. 6) to deter- 
mine whether or not substantial registration between the 
working set image sets exists. Where substantial regis- 
tration does exist, control passes out of the process and 
the process ends. However, at block 156, where sub- 

55 stantial registration between the working set image sets 
does not exist, control passes to block 158 where com- 
puter 60 identifies first and second adjacent image sets 
within the working set that are misregistered to the 
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greatest degrees. For instance, in the present example, 
computer 60 may identify working set image sets 110 
(v N ) and 112(v N ) as the most misaligned or misregis- 
tered adjacent image sets. Herein it will be assumed that 
image set 110(v N ) and image set 112(v N ) are identified 
as the first and second sets, respectively. 
[0054] Continuing, at block 1 60 computer 60 identifies 
the phase image set corresponding to the same diastolic 
phase as the second image set that most closely regis- 
ters with the first image set as a substitute image set. 
Referring once again to Fig. 8, in the present example, 
computer 60 identifies phase image set 1 1 2(v N+1 ) as the 
image set corresponding to the second phase image set 
1 1 2 (v N ) that is most closely aligned with the first image 
set 110(v N ). At block 162, computer 60 substitutes the 
substitute image set for the second image set generat- 
ing a modified working set. Thereafter, control again 
passes back up to block 156 where computer 60 con- 
tinues to determine whether or not the working set im- 
age sets are substantially registered and, where they 
are not substantially registered, continues to substitute 
other phase images sets until substantial or optimal reg- 
istration results. 

[0055] It should be understood that the methods and 
apparatuses described above are only exemplary and 
do not limit the scope of the invention, and that various 
modifications could be made by those skilled in the art 
that would fall under the scope of the invention. 
[0056] For example, there are many different image 
comparison algorithms that measure the degree of mis- 
registrationin the prior art that could be used to perform 
the image set comparisons and the present invention is 
not limited to any one such process. 
[0057] As another example, while the present inven- 
tion as described above is one wherein raw helical data 
is stored and image sets are only generated after a sys- 
tem operator selects a location SP (see again Fig. 4) for 
image processing, it should be appreciated that other 
embodiments exist where, for example, phase depend- 
ent two-dimensional image sets may be generated for 
each diastolic phase and the phase dependent image 
sets may be stored prior to selection of location SP. In 
this case, after a location SP is selected, computer 60 
may simply access the phase image sets for each di- 
astolic phase and select the image sets from each di- 
astolic phase that are most closely aligned with location 
SP. Thereafter, the process described in Fig. 8 including 
blocks 156, 158, 160 and 162 could be performed as 
described above. 

[0058] In addition, while the inter- diastolic phase 
ranges are shown above as overlapping, the present in- 
vention may be performed with truly sequential phase 
ranges where adjacent phase ranges do not overlap or, 
in the alternative, may be performed with phase ranges 
that are separated by some small range. Moreover, 
while the invention is particularly useful for image pur- 
poses during diastolic phases, clearly, other phases 
such as the systolic phases may be targeted for image 



processing according to the present invention. In fact, 
because the systolic phases is generally shorter than 
the diastolic phase the present invention may be partic- 
ularly advantageous for imaging during the systolic 
5 phase. Furthermore, while the invention as described 
above in the context of EKG gated applications, clearly 
the present invention is useful other applications where 
periodic structural motion occurs. 

[0059] For the sake of good order, various aspects of 
10 the invention are set out in the following clauses:- 

1 . A method to be used with a CT imaging system 
for minimizing phase misregistration among the two 
dimensional CT images that make up a working im- 

15 age set corresponding to a region of interest (ROI) 
(120), the method comprising the steps of: 

identifying (148) separate sequential acquisi- 
tion phases during an acquisition period; 
20 obtaining (142) helical data corresponding to 

the separate acquisition phase; 
for each acquisition phase: 

(i) identifying (148) at least two different 
25 phase ranges; and 

(ii) generating (152) and storing images as 
phase image sets corresponding to the dif- 
ferent phase ranges; 

30 comparing (156) the phase image sets corre- 

sponding to sequential ranges to identify sub- 
stantially aligned phase image sets; and 
selecting (158, 160, 162) the substantially 
aligned phase image sets as the working image 
35 set for the ROI. 

2. The method of clause 1 wherein the step of iden- 
tifying acquisition phases includes the step of mon- 
itoring the cycle of a cyclical anatomical function 

^0 and wherein the step of obtaining includes gating 
helical data acquisition so as to coincide with the 
cycles of the anatomical function. 

3. The method of clause 2 wherein the anatomical 
^5 function is a heart beating cycle and each acquisi- 
tion phase corresponds to at least one heart cycle 
phase. 

4. The method of clause 3 wherein the at least one 
50 heart cycle phase is the diastolic phase. 

5. The method of clause 3 wherein the at least one 
heart cycle phase is the systolic phase. 

55 6. The method of clause 1 wherein the step of ob- 
taining includes the steps of acquiring (144) a full 
set of helical data while monitoring the cycle of a 
cyclical anatomical function and storing the helical 
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data as a function of the cyclical anatomical 
the step of identifying acquisition phases includes, 
after the helical data is stored, identifying a separate 
phase (146) of the anatomical function cycle and 
the step of obtaining further includes the step of re- 5 
trieving the helical data corresponding to the sepa- 
rate acquisition phase. 

7. The method of clause 6 wherein the anatomical 
function is a heart beating cycle and each acquisi- 10 
tion phase corresponds to at least one heart cycle 
phase. 

8. The method of clause 1 further including the step 

of, prior to comparing, identifying a specific phase '5 
point from within an acquisition period and wherein 
the step of comparing includes selecting phase im- 
age sets most closely corresponding to the identi- 
fied phase point as a working set and comparing 
the working set image sets to identify misalignment 20 
and wherein the step of selecting includes substi- 
tuting non-working set image sets for working set 
image sets where the substitutions reduce mis- 
alignment and provide a modified working image 
set. 25 

9. The method of clause 8 wherein the step of com- 
paring the working set images includes the steps of: 

(a) identifying (1 58) working set image sets cor- so 
responding to sequential acquisition phases 
that have the greatest degree of misalignment 

as first and second image sets; 

(b) identifying (160) a non-working set image 

set corresponding to the same acquisition 35 
phase as the second image set as a substitute 
image set; 

(c) substituting (162) the substitute image set 
for the second image set; and 

(d) repeating steps (a) through (c) until a sub- 40 
stantially aligned working image set results. 

1 0. The method If clause 9 wherein the step of iden- 
tifying a non-working set image set includes the 
step of identifying a non-working set image set cor- 45 
responding to the same acquisition phase as the 
second image set that is most aligned with the first 
image set as the substitute image set. 

11 . The method of clause 1 wherein the phase rang- so 
es overlap. 

12. The method of clause 1 wherein the phase rang- 
es are consecutive. 

55 

13. The method of clause 1 wherein the step of iden- 
tifying at least two phase ranges includes the step 
of identifying more than two phase ranges for each 



1 4. An apparatus to be used with a CT imaging sys- 
tem for minimizing phase misregistration among the 
two dimensional CT images that make up a working 
image set corresponding to a region of interest 
(ROI), the apparatus comprising: 

a processor (60) running a pulse sequencing 
program to perform the steps of: 

identifying (148) separate sequential ac- 
quisition phases during an acquisition pe- 
riod; 

obtaining (142) helical data corresponding 
to the separate acquisition phase; 
for each acquisition phase: 

(i) identifying (148) at least two differ- 
ent phase ranges; and 

(ii) generating (152) and storing imag- 
es as phase image sets corresponding 
to the different phase ranges; 

comparing (156) the phase image sets cor- 
responding to sequential ranges to identify 
substantially aligned phase image sets; 
and 

selecting (158, 160, 162) the substantially 
aligned phase image sets as the working 
image set for the ROI. 

15. The apparatus of clause 14 wherein the pro- 
gram causes the processor to perform the step of 
identifying acquisition phases by monitoring the cy- 
cle of a cyclical anatomical function and wherein the 
step of obtaining includes gating helical data acqui- 
sition so as to coincide with the cycles of the ana- 
tomical function. 

16. The apparatus of clause 15 wherein the ana- 
tomical function is a heart beating cycle and each 
acquisition phase corresponds to at least one heart 
cycle phase. 

1 7. The apparatus of clause 1 (1 4) wherein the pro- 
gram causes the processor to perform the step of 
obtaining by acquiring (1 44) a full set of helical data 
while monitoring the cycle of a cyclical anatomical 
function and storing the helical data as a function of 
the cyclical anatomical data, the program causes 
the processor to perform the step cf identifying ac- 
quisition phases by, after the helical data is stored, 
identifying a separate phase (1 46) of the anatomical 
function cycle and wherein the the step of obtaining 
further includes retrieving the helical data corre- 
sponding to the separate acquisition phase. 
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18. The apparatus of clause 14 wherein the pro- 
gram further causes the processor to, prior to com- 
paring, identify a specific phase point from within an 
acquisition period and wherein the program causes 
the processor to perform the step of comparing by 
acquisition phases by selecting phase image sets 
most closely corresponding to the identified phase 
point as a working set and comparing the working 
set image sets to identify misalignment and wherein 
the program causes the processor to perform the 
step of selecting by substituting non-working set im- 
age sets for working set image sets where the sub- 
stitutions reduce misalignment and provide a mod- 
ified working image set. 

19. The apparatus of clause 14 wherein the phase 
ranges overlap. 

20. The apparatus of clause 14 wherein the phase 
ranges are consecutive. 

21 . The apparatus of clause 14 wherein the step of 
identifying at least two phase ranges includes the 
step of identifying more than two phase ranges for 
each acquisition phase. 



Claims 

1 . A method to be used with a CT imaging system for 
minimizing phase misregistration among the two di- 
mensional CT images that make up a working im- 
age set corresponding to a region of interest (ROI) 
(120), the method comprising the steps of: 

identifying (148) separate sequential acquisi- 
tion phases during an acquisition period; 
obtaining (142) helical data corresponding to 
the separate acquisition phase; 
for each acquisition phase: 

(i) identifying (148) at least two different 
phase ranges; and 

(ii) generating (152) and storing images as 
phase image sets corresponding to the dif- 
ferent phase ranges; 

comparing (1 56) the phase image sets corre- 
sponding to sequential ranges to identify sub- 
stantially aligned phase image sets; and 
selecting (158, 160, 162) the substantially 
aligned phase image sets as the working image 
set for the ROI. 

2. The method of claim 1 wherein the step of identify- 
ing acquisition phases includes the step of monitor- 
ing the cycle of a cyclical anatomical function and 
wherein the step of obtaining includes gating helical 



data acquisition so as to coincide with the cycles of 
the anatomical function. 

3. The method of claim 1 or 2 wherein the anatomical 
5 function is a heart beating cycle and each acquisi- 
tion phase corresponds to at least one heart cycle 
phase. 

4. The method of claim 3 wherein the at least one heart 
10 cycle phase is the diastolic phase. 

5. The method of claim 3 wherein the at least one heart 
cycle phase is the systolic phase. 

15 6. An apparatus to be used with a CT imaging system 
for minimizing phase misregistration among the two 
dimensional CT images that make up a working im- 
age set corresponding to a region of interest (ROI), 
the apparatus comprising: 

20 

a processor (60) running a pulse sequencing 
program to perform the steps of: 

identifying (148) separate sequential ac- 
25 quisition phases during an acquisition pe- 

riod; 

obtaining (142) helical data corresponding 
to the separate acquisition phase; 
for each acquisition phase: 

30 

(i) identifying (148) at least two differ- 
ent phase ranges; and 

(ii) generating (1 52) and storing imag- 
es as phase image sets corresponding 

35 to the different phase ranges; 

comparing (156) the phase image sets cor- 
responding to sequential ranges to identify 
substantially aligned phase image sets; 
40 and 

selecting (158, 160, 162) the substantially 
aligned phase image sets as the working 
image set for the ROI. 

45 7. The apparatus of claim 6 wherein the program 
causes the processor to perform the step of identi- 
fying acquisition phases by monitoring the cycle of 
a cyclical anatomical function and wherein the step 
of obtaining includes gating helical data acquisition 
so so as to coincide with the cycles of the anatomical 
function. 

8. The apparatus of claim 6 or 7 wherein the anatom- 
ical function is a heart beating cycle and each ac- 

55 quisition phase corresponds to at least one heart 
cycle phase. 

9. The apparatus of claim 6, 7 or 8 wherein the pro- 
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gram causes the processor to perform the step of 
obtaining by acquiring (144) a full set of helical data 
while monitoring the cycle of a cyclical anatomical 
function and storing the helical data as a function of 
the cyclical anatomical data, the program causes 5 
the processor to perform the step of identifying ac- 
quisition phases by, after the helical data is stored, 
identifying a separate phase (146) of the anatomical 
function cycle and wherein the the step of obtaining 
further includes retrieving the helical data corre- io 
sponding to the separate acquisition phase. 

10. The apparatus of any one of claims 6 to 9 wherein 
the program further causes the processor to, prior 
to comparing, identify a specific phase point from *5 
within an acquisition period and wherein the pro- 
gram causes the processor to perform the step of 
comparing by acquisition phases by selecting 
phase image sets most closely corresponding to the 
identified phase point as a working set and compar- 20 
ing the working set image sets to identify misalign- 
ment and wherein the program causes the proces- 
sor to perform the step of selecting by substituting 
non-working set image sets for working set image 
sets where the substitutions reduce misalignment 25 
and provide a modified working image set. 
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